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NUMERICAL EXPLORATION OF A HEXAGONAL STRING 

BILLIARD 

HANS L. FETTER 



Abstract. In this paper we are interested in the motion of a ball inside a 
billiard table bounded by a particular smooth curve. This table belongs to a 
family of billiards which can all be drawn by a common process: the so-called 
gardener's string construction. The classical elliptical billiard is, of course, the 
foremost member of this family. So it should come as no surprise that our 
hexagonal string billiard shares many basic properties with the latter, but, on 
the other hand, also exhibits some essential differences with it. 



(N 

a,; 

j~l 1. Background 

-(— > \ 

Let us consider the motion of a point inside a plane billiard table bounded by a 
closed convex curve. This point will always move along a straight line until it hits 
the boundary where it is reflected according to the well known principle: the angle 
of reflection is equal to the angle of incidence. 
^. \ In order to get a billiard table with a sufficiently smooth boundary we can use 

<3\ ■ the familiar method known as a the gardener's or the string construction. First 

we need to choose some convex polygon K and then we proceed to wrap a loop of 
inelastic string around it, pulling the string tight at a point P and then moving this 
point P around. Note that when this technique is applied to a closed line segment 
one obtains an ellipse. Much is known about the billiard problem inside an ellipse: 



en 

~"» \ see for instance the books by Chernov and Markarian[6] and Tabachnikov 35 , and 

also the articles by Berry 4 , Korsch and ZimmerpT] and AcquistapacefT] just to 
name a few. 

Another choice for the convex polygon which has received a fair amount of at- 
tention is that when K is an equilateral triangle. References include Hubacher |18j . 



Gutkin and Katok [13], Gutkin and KnillpT], Turner [36]. Here again the result is 
a convex domain whose boundary is comprised of portions of ellipses. 

In Fig. [T]we illustrate the gardener's construction of a table $7 for the case when 
if is a regular hexagon. The point P, as it moves around, traces out a curve dfl 
which again consists of a certain number of elliptical arcs. 

Varying the length I of the inelastic string one gets, of course, a whole family 
of billiard tables. We need to restrict our attention, however, to only one special 
value of I. This choice is dictated by the smoothness we get for the boundary <9f2. 
So it is for it only for this particular case that we want to obtain a description, 
as complete as possible, of all the trajectories as the billiard ball (point particle) 
bounces elastically off the walls. 

In closing this introduction let us just mention that the study of billiards and in 
particular the subject of billiards inside smooth convex curves is a very active area 
of research in the field of dynamical systems. The interested reader may consult 

l 
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Figure 1 . Hexagonal string billiard 



some of the following books and research articles dealing with various important 
aspects of this topic (see [2D1 [221 El ES] and [H [T31 HH [HI HH H31 HU respectively). 



2. A FAMILY OF BILLIARDS OF CLASS C 2 

The billiard table obtained by the string construction for an equilateral triangle, 
regardless of the string length used, is of class C 1 (see Hubacher [IB])- I* 1 general, 
this is also true for all the other billiard tables obtained from regular polygons by 
the string construction, including the regular hexagon. However, with a proper 
choice of the string length I one can construct a table having greater regularity. 
To accomplish this let us first outline the setup and also the strategy behind the 
construction process for a whole family of "smooth" string billiards. To start we 
select K as any convex regular n-gon (n > 5). So right from the beginning we 
exclude both the equilateral triangle and the square from further considerations. 
The reason for that is because we need to be able to stellate K. A stellated polygon 
can be derived from a regular polygon by adding identical (congruent) isosceles 
triangles to all its sides. Instead of adding them all we add only one of those 
triangles to our given polygon K. The perimeter of this resulting polygon is the 
sought for special value for the string length I. Shortly we shall find an explicit 
expression for it. A portion of a regular polygon having vertices F\ , F% , . . . , F n with 
side length 2 is shown in left part of Fig.[2j The coordinates of the vertices, in terms 
of the interior angle a = — -tr, are as indicated below: 
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Fi = [-1,0] 
Fa = [1,0] 

F3 = [1 — 2 cos a, —2 sin a] 

F4 = [—1 — 2 cos a + 4 cos 2 a, 2(— 1 + 2 cos a) sin a] 

F n _x = [1 + 2 cos a — 4 cos 2 a, 2(— 1 + 2 cos a) sin a] 
F n = [— 1 + 2 cos a, —2 sin a] 

Now let us add the isosceles triangle with top vertex G\ to the side Fi,F% 
(see left part of Fig. [5]). In this figure we have actually added several identical 
(congruent) isosceles triangles to the sides of the regular n-gon. They have vertices 
G\, G2, • • • j G„, whose coordinates are given below: 



G 2 



0, 



sin a 



cos a 
cos a — 1 
cos a 



Un — 



cos a — 1 



.0 



cos a 



Each of the equal sides of the isosceles triangle F1F2G1 has length d = |i*iGi| = 
|F 2 Gi| =-l/cosa. 

So finally we get the expression for the particular string length we were looking 
for: 

I = 2(n - 1) + 2d = 2(n - 1) rlzr- 

cos(^7r) 

So much for the setup. In the right part of Fig. [2] we illustrate the result of the 
construction process described above for the case when K is a regular nine-gon. 
That we get billiard tables with greater regularity than just G 1 is the contents of 
the following theorem. 

Theorem 2.1. The generalized gardener's construction with length of the string 

2 
Z = 2(n-1)- 



cos (^7r) 

provides us with a family of billiard tables (one for each convex regular n-gon, n > 5) 
which are globally G 2 . 

For the proof see O 

Some immediate consequences from this construction process include the follow- 



ing: 



As opposed to the general case we need to concern ourselves with only one 
kind of ellipse. 

The closed boundary curve we get is a smooth union of n elliptical arcs. In 
the general case one gets twice that number. 
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Figure 2. Construction of n-gonal string billiard and enneagonal 
string billiard 



• Since an ellipse will always have strictly positive curvature the curvature k 
is a continuous, strictly positive function on d£l. 

From this point on we shall discuss one specific billiard only: the hexagonal 
string billiard. We shall now provide some specific details for it (see Fig. [3]). Let K 
be the regular hexagon with sidelength 2. Its vertices have the coordinates 



F 4 



(-1,V3), F 2 



(1,V3), 



^3 

F 6 



(2,0) 
(-2,0) 



Since n = 6, a = -£ and d 



-1/ cos a = 2 the length of the string I — 2(n — 1) + 
2d = 14. Each of the six elliptical arcs uses a unique pair of foci Fi, Fj , where i and 
j are either both even or both odd. To reflect this in our notation, we have used 

arc ij for an arc focused at Fi , Fj . 

The equation for the arc 24 is 

(z-i) 2 , V 2 



G 



9 



1 



whereas the equation for the arc 13 is given by 

lly 2 + 2v / 3(a: - 6)y + 9x 2 - Ylx = 60 

At the common point (3, y/3) the first derivative for both yields — \/3 while the 
second derivative is — |Vo- 

As for the curvature k of any elliptical arc, a simple calculation gives: 

V6 ^ ^ 3V3 
9 - " 16 
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M)\(3,V3) 




cos (7) = 



r + 9 

27 -f 



(s-1) 2 , y 2 = 1 

6 9 



^2(27-y 2 ) 3/2 



54 



Figure 3. Data of hexagonal billiard 



3. Types of orbits 

The orbits of billiards are sequences of line segments on the billiard table with 
the segments corresponding to the straight paths which the ball must follow within 
the region until it hits the boundary where it reflects elastically. Each such orbit is 
determined by the choice of the initial location for the billiard ball as well as the 
direction for the shot. 

It is well-known that in a generic or typical billiard one can identify different 
types of trajectories: periodic, quasi-periodic (librational), whispering-gallery (ro- 
tational) and chaotic. We shall try to describe and illustrate them in the context 
of the hexagonal string billiard. 



3.1. Periodic orbits. In almost any study on billiards initial interest centers on 
periodic orbits: those in which the particle follows the same path endlessly. In 
general these paths can be recognized easily because they simply correspond to 
inscribed simple polygons or star polygons. When, for example, a convex billiard 
table has two perpendicular symmetry axes then we immediately get a period four 
orbit (Hasselblatt and Katok [E]). Our billiard has six axes of symmetry, so we 
get twelve points from their intersection with the boundary of the billiard table. 
Selecting arbitrarily any of these points as the initial point and then connecting it 
to the next one by skipping over exactly fc — 1 = 0, . . . , 5 intermediate points we 
are bound to get a closed figure. Since this figure is always invariant under the 
symmetries, it follows that at each point of intersection with the table the angle 
of incidence is equal to the angle of reflection. Now, according as k = 1,...,6 
the corresponding periodic orbits will close onto themselves after 



12 

(12, k) 



bounces, 
respectively, where (n, k) denotes the greatest common divisor of n and k. Some of 
these basic periodic orbits for our billiard are shown in Fig. 01 

When dealing with periodic trajectories a very natural question is: how many 
are there performing n bounces and at the same time going around the boundary 
k times? Birkhoff [5] provided an astonishing answer to the previous question as 
early as 1927 (see also Tabachnikov [35]): 
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(a) Unstable orbits of (b) Stable orbits of pe- (c) Unstable orbits of pe- 
period 2, 3, and 6 riod 2, 3, and 6 riod 4, 12, and 12 



Figure 4. Periodic orbits j 1 ^} where k = 1 . . . 6 



Theorem 3.1. For any n > 2 and every fc < §, relatively prime, there exist two 
geometrically distinct n— periodic trajectories with the rotation number k. 

So, for example, the 2-, 3- and 6-bounce orbits in Fig. HJa) and those in Fig.HJb) 
are geometrically distinct. They are similar though, with a similarity ratio of 



2 v / 3(-l + \/6) 



1.0042 



From BirkhofFs theorem we are lead to the conclusion that there are always 
at least <p(n) distinct n— periodic trajectories, where <p(n) denotes Euler's totient 
function. Thus when n = 17 we must have ^(17) = 16 distinct 17— periodic tra- 
jectories. The reader can appreciate the eight unstable 17-sidcd star polygons that 
we constructed for the hexagonal string billiard in Fig. [SJ 




Figure 5. Star polygons: {\ } ■ ■ ■ { J} 
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It is clear that, because our billiard has six- fold symmetry, we can rotate these 
geometrically distinct n— periodic trajectories by certain angles. The order of rota- 
tional symmetry is given by , e 6 n s and so this yields a grand total of 

6 



(6,n) 
n— periodic trajectories for each n > 3. 



tp(ri) 



The first few values of this expression are given in the following Table Q] 



n 


Teiky^W 


n 


Tg^yvC") 


3 


4 


8 


12 


4 


6 


9 


12 


5 


24 


10 


12 


6 


2 


11 


60 


7 


36 


12 


4 



Table 1. Total number of n— gons. 



3.2. Quasi-periodic orbits. Sometimes one also encounters certain orbits which 
explore only restricted parts (segments) of the boundary <9f2. Every (stable) simple 
polygon or star polygon (the two-bounce orbit included) gives rise to an infinite 
number of these quasi-periodic orbits. They resemble the underlying polygons: it 
is as if these just fattened up. So that's why we will also refer to them simply as fat 
(star) polygons. Orbits which explore 2,3,4,6 and 9 parts can be seen in Figs. [5] 
and [7] More examples can be found in Figs. [TU1 H4l 1221 and l24l 




Figure 6. A typical bouncing orbit or fat polygon of type {^} 



3.3. Whispering gallery orbit. Another kind of orbits, which on occasions are 
also referred to as skipping trajectories (see Berry 0], Chernov and Markarian[B]), 
are very common for most billiards. They are characterized by the fact that they 
bounce all round d£l densely filling a ring-like shaped region. 

There are many of these orbits in our billiard and we show several of them in 
Fig. [8] 
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Figure 7. Quasi-periodic orbits or fat polygons of type: 
{.MM} -d {»} 




Figure 8. Whispering gallery type orbits 



3.4. Chaotic orbit. Most billiards exhibit chaotic behavior to some extent (circles 
and ellipses excluded). Chaotic trajectories look disordered, erratic and there is no 
discernible pattern. They usually manifest themselves more clearly in the surface 
of section, which we will introduce later. For the moment let us just say that there 
is no evidence that our billiard has these kinds of orbits. 
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4. Some analogous properties of the elliptical and the hexagonal 

string billiard 

One of the most important properties of elliptic billiards is the following result 
which essentially tells us that we can group the trajectories into three families. 
A trajectory for the billiard inside an ellipse either 

(1) always passes through the two foci alternately (focal trajectory), or 

(2) always intersects the open segment between the two foci (inner trajectory), 
or 

(3) never intersects the closed segment between the foci (outer trajectory). 

(See Siburg [32], Tabachnikov [35], Chernov and Markarian [6]). 
In a completely analogous fashion for the hexagonal string billiard we also have 
three families of trajectories: focal, inner, and outer. 

Theorem 4.1. Consider a trajectory issuing from the point Pq inside a hexagonal 
string billiard with vertices Pi , P 2 % ■ ■ ■ ■ 

i) If the segment PqP\ of the billiard trajectory is on a supporting line of K , 
then every segment PiPi + \ will be on a supporting line of K . 

ii) // the segment PqPi of a billiard trajectory does intersect the hexagon K 
then all segments PiPi + i will intersect it. 

hi) // the segment PqP\ of a billiard trajectory does not intersect the hexagon 
K then no segment PiPi+i will intersect it. 

A sketch for the proof of this result can be found in [B] 

Another important result for the billiard inside an ellipse states that a billiard 
trajectory through the foci converges to the major axis of the ellipse. (See for 
example Tabachnikov [35], Batschelet [3], Frantz [10], Hasselblatt and Katok [16] . 
Wilker [37], Moser and Zehnder [27]). 

In practice this means that any focal orbit quickly becomes indistinguishable 
from a repeated tracing of the major axis of the ellipse. 

Whereas a focal orbit in the case of the ellipse converges to the "to and fro" orbit 
(the unstable two-bounce orbit) , in that of the hexagonal string billiard it converges 
to either one of the equilateral triangles (the unstable three-bounce orbits). 

Before stating the result in more precise terms we need an auxiliary result which 
provides bounds for the sizes of the angles between incoming and outgoing rays in 
focal orbits. (See Fig. |3j. 

Theorem 4.2. For a focal orbit of the hexagonal string billiard the angle between 
incoming and outgoing segments is limited to the range 

60° < 7 < 70.53° 

Proof in [C] 

Without loss of generality let us choose the three-bounce orbit shown in Fig. S(a) . 
The vertices for this triangle are £1 = (0, 2-\/3), £3 = (3, — \/3), and S5 = (—3, — \/3). 
Consider a billiard trajectory, passing through the focus F2, whose initial point Pq 

is on arc 24. See Fig. [9] Observe that, for example, point P4 is closer to point Si 
than Pi, and then P7 is even closer to it and so on. Already we hnd that PrPsPg 
resembles triangle S1S5S3 very much. 
Under these circumstances we have: 
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Theorem 4.3. If a billiard trajectory travels along a supporting line of K , so that 
it goes alternately through the focal points F<z,F6,F& then its trajectory will quickly 
become indistinguishable from a repeated tracing of the triangle S1S5S3. 



See [Dj for proof. 
















5tC 

/\y 

/ 1/ 

i / 


T7\~ 


\s 


\ 
\ 


/ s J 




T/ 1 


\ 


»T1 






\\ 1 


/ L 

i Ir 






\ S 3 1° 


i A 






\\ ' 


I f\ 






J 

\ ■ 
\ ■ 


\J\ 




" F 4 


/ 9 


S K 


-<" 


3> 


/ 


^ — 




^,0 


y 


*\ 






^ 



(a) Trajectory approaching the triangle 
S1S5S3 



( b ) {*i}>{¥>i}. {«;} 



Figure 9. Focal orbits 

5. Stability 

Several authors include in their work a brief discussion on the stability of certain 
n— periodic orbits through the use of the so called deviation matrix M. We also 
want to mention something about stability. For this purpose we have chosen the 
j 1 ^} orbit, which can be seen on the top left part of Fig. [TO] This orbit is not the 
same as the earlier one we came across in Fig.@|c) (unstable j 1 ^} orbit). 

So following Berry [4], Lopac et al. [25], [21], and Robnik [3T] we have to 
construct the deviation matrix M corresponding to the selected orbit: 

M = Mi )2 M 2 , 3 • • • M 12 ,i 
where each 2x2 matrix M$ & can be expressed as 



Pik 



M ilk = 



RiRk 



Ri 



j Sill Qfc 

+ 



pik 



sin a 
sin_Qfc 
sin Oii 



+ 



sin afc 
Pik 



i?,^ sin oti 



The symbol Ri (Rk) denotes the radius of curvature at the point of impact Pi 
(Pk), on (offe) the angle which the departing ray makes with the tangent on the 
billiard boundary at the point of impact Pi (Pk) and pik the length of the path 
between two consecutive impact points Pi and Pk- 

For our particular choice we only get two kinds of these 2x2 matrices M^ since 
only the path lengths change between impacts, all other data remaining constant. 

In fact, if we define 



3a/(4947-2328V3) 



97 
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then we get 

(1) R l = R k = ^(w^t3^' 2 

(2) a.i — ak = 5-7r/12, so that sin(ai) = sin(afc) = — ^ 

(3) the two values which alternate 



Plk = J(fr 2 + 8V^r 2 + 8r + 6\/3t + 3) 



8-r-L. 4V3, 



So 



r 



Mi, a = M 3 , 4 = M 5;6 = M 7 , 8 = Af 9 ,io = M u ,i 2 = T 
M 2 , 3 = M 4 , 5 = M 6 , 7 = M 8)9 = A/l0,ll = Mi2,l = S 

0.9830565623814575557 -7.1795378524870580504 
0.00467993843741796985 0.9830565623814575557 

0.9700724323780286883 -7.1325295852446540944 
0.00826627849706319570 0.9700724323780286883 
Therefore our deviation matrix M turns out to be 

-.87488446006857883011 -17.010433690097981202 



where 



and 



U irSl l .015415741802173798979 -.84327883276452294763 

Orbital stability depends on the eigenvalues of M. In fact, if by tr(M) we denote 
the trace of M then it is known that (Berry [3]) we have the following possibilities: 

(1) if |tr(M)| < 2 then the orbit is stable 

(2) if |tr(M)j > 2 then the orbit is unstable 

(3) if |tr(M)| = 2 then the orbit is neutrally stable. 
For our case we simply get that 

tr(Af) = -1.7181632928331017777 

which guarantees the stability of our chosen orbit. This analytically verified sta- 
bility can also be observed numerically: see Fig.[TUl where in addition to the { 12 } 
polygon we show several quasi-periodic orbits. 

6. Forbidden inner region 

In many of the figures (see for example Figs. fTO)) we can discern a closed 
convex proper subset J- of tt which lacks segments of the corresponding billiard 
trajectory. We shall call it a "forbidden inner region", a term which we have bor- 
rowed from Korsch and Zimmer 21j. It is separated from the allowed region by a 
non diffcrcntiable piecewise smooth closed curve. By all of this we are, of course, 
reminded of the notion of a caustic. There are several different meanings of the 
word caustic in billiards. The one that seems the most appropriate for us is (see 
Gruber Q2J, Klee and Wagon [15], Turner [3B]): 

A convex caustic for a convex table O is a closed convex proper subset C of fi 
such that whenever the initial segment of a billiard path in il lies in a supporting 
line of C, then every later segment also lies in a supporting line of C. All billiard 
tables obtained by the string construction have a special caustic. The hexagon K 
is such a caustic and we already established this in Theoreir l4.ll 
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Figure 10. Stable 



{ 5 ) orbit and several fat { 5 } 's 



Are the forbidden inner regions we get in the other cases caustics as well? For 
us to have a convex caustic every line segment of a trajectory needs to touch 
C. However careful inspection of the billiard trajectories reveals that a few line- 
segments do not touch C (see Fig. ITT]) . 

We need to introduce an alternative description for the forbidden inner region. 
We define the forbidden inner region T as the intersection of all the "left" half- 
planes to the lines containing the segments of a billiard trajectory (see Hasselblatt 
and Katok [To]). 

We have included Fig. [TTJ which allows us to follow the construction process of 
a forbidden inner region with an increasing number of half-planes. 

A problem remaining is that of trying to find an analytical expression for the 
piecewise smooth closed curve bounding the forbidden inner region. The orbits 
shown in Figs. [7] [5] Q7J and [12] strongly suggest that a hexagon with curved sides 
might be appropriate for the curve we are looking for. 

Using elliptical arcs for the sides of the curved hexagon we get the results shown 
in Fig. [T3] and also in Fig. [TTJ 

7. POINCARE SECTION 

So far we have presented all orbits as merely sequences of directed line segments 
in the configuration space. This approach, even if we observe the evolution of a 
large number of these trajectories, gives us only a limited understanding of the 
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Figure 1 1 . Construction of forbidden inner region 




Figure 12. Forbidden inner regions for groups of quasi-periodic orbits: 
6 {D's, 6 {D's, 2 {*}',, and 3 {*}'s 
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Figure 13. Arcs of ellipses separating forbidden inner region 




Figure 14. Piecewise smooth curve separating forbidden inner region 
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complicated dynamics associated with the motion. A more convenient way to visu- 
alize the dynamics is through the use of Poincare's surface of section (SOS), which 
enables us to display the character of each particular trajectory. The surface of 
section, see for instance Sussman and Wisdom[34j, can be described by means of 
two coordinates in the 

(1) horizontal direction: arc length s from an arbitrarily chosen fixed point O 
on the boundary of the billiard, and 

(2) vertical direction: the angle 9 between positive tangent direction and out- 
going or reflected billiard ray at a point of impact. 

This then defines the following Poincare section 

V:={(s,e)\se[0,\dn\],9e[0,n]} 

There are, qualitatively speaking, three different classes of trajectories possible 
on a SOS, clearly differentiated by the dimension d of the subspace of the section 
that they explore (Berry [I] .HavlifTT] .Korsch and Jodl[2"0] .Sussman and Wisdom[3~4"]). 
They might 

(1) generate a finite set of k discrete points corresponding to a periodic orbit 
having period k (d = 0) 

(2) eventually fill out an invariant curve, which may either run from edge s = 
to edge s — 2ir as an undulating line corresponding to a whispering-gallery 
or a focal orbit, or consist of a group of ovals or islands corresponding to a 
quasi-periodic orbit (d = 1) 

(3) scatter over a region eventually filling a whole area corresponding to an 
irregular or chaotic trajectory (d = 2). 

Because our billiard has the same symmetry as a hexagon, we do not need to 
plot the entire surface of section. We shall usually only display the upper half of it 
and note the fact that the SOS is periodic with period six. 

To give the reader an idea of the structure of the SOS we first present an image 
in which we have included points and curves corresponding to some 70 orbits, each 
followed through 480 bounces (see Fig. [T5j) . In Fig. [16] we do not show all of them 
at once as was done before, but rather in small groups, for different ranges of the 
angle 9. 

Several close-ups of the six-bounce orbits which can clearly be recognized in the 
plot on the right in the top row of Fig. [Jj)] are shown in Fig. 1171 

Another close-up (Fig. ITS)) , this time of the bottom plot from Fig. [TBI shows both 
undulating lines and ovals. The big ovals surround the invariant points correspond- 
ing to the stable periodic orbits associated with { 5 } 3 {23} ' {3}' an< ^ {13}- 

Finally, we decided also to include the curve generated in the surface of section 
by a large number of focal orbits (see Fig.[l9|). It is interesting to note that in this 
case it is even possible to get an explicit expression for each arc that appears. The 
coordinates of the points on a reference arc for the focal curve are given by 



where 



E f y/(l + 3cos(2y))(-l + cos(2y)) V3\ 
{ l-oos(2y) ' 3 J ,y , 



arccos(— 1/3) arccos(— 1/2) 

< y < 



2 ~" - 2 

and E(x, k) is simply an elliptic integral of the second kind (see Davis [?])• 
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Figure 15. Surface of section: 70+ orbits followed through 480 bounces 



8. Additional periodic orbits and some comments on the structure 
OF ORBITS in configuration space 

All the periodic orbits we have presented so far are those corresponding to star 
polygons of the kind {£}, with 1 < k < n/2 and (n, k) = 1, which always come in 
pairs: one stable and the other one unstable. Moreover these orbits are isolated, as 
is customary for a generic billiard (Berry [4]). 

So, now, a natural question arises as to whether there are any other periodic 
orbits besides those just mentioned. After all, many authors, among them Berry [4], 
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Figure 16. Surface of section: same as Fig. [15j but more detailed 



Korsch and Zimmer[2"Tj, Okai et al. [28] . and Sieber (33J include periodic orbits in 
their work which do not correspond to star polygons. Of course we found many 
additional periodic orbits for our hexagonal string billiard of a different type, some 
of which are shown in Figs. [20] and [21] 

In the captions of these figures we have used the notation { ^} in spite of the fact 
that n and k are no longer relatively prime, because we simply want to express the 
basic fact that the orbit has period n and that it winds around the table k times. 
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Figure 17. Surface of section: close-up of {^} orbit 



Interestingly enough, all of these orbits are again isolated, but, in contrast with 
the ones corresponding to star polygons, these additional periodic orbits turn out 
to be neutrally stable! 

Now let us include a brief comment about the organization of periodic and quasi- 
periodic trajectories in configuration space. 

Periodic trajectories inside an elliptic billiard are very well organized. Recall 
that for fixed n and fc, relatively prime positive integers with k < ro/2, there 
is a continuous family of periodic orbits of the form {^}. All the trajectories 
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Figure 18. Surface of section including fat j 1 ^}, {23} , {3}, and{'J°} 



corresponding to such a family are arranged around a common ellipse or hyperbola 
(see Poncelet's theorem [35]). 

All periodic orbits for our billiard are isolated. So for fixed n and k, where n 
and k are relatively prime positive integers with k < n/2, we want to consider 
not only the periodic orbits of the form |?| but all the quasi-periodic orbits or fat 
polygons {^} as well. If there are any periodic orbits of the form { m }, j < m/2 with 
m/j = n/k we need to include them too. We want to see how they are organized in 
configuration space: if they are arranged around a common forbidden inner region. 

To begin we will show what "goes on between two consecutive" unstable periodic 
orbits: {{} and {J}. 

The claim now is that for either n = 3 or 4 the pair of {"} unstable periodic orbits 
"encloses" a {"} stable periodic orbit, periodic orbits of the form { m }, j < m/2 
with m/j = n and fat polygons {"} as well. 
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Figure 19. Surface of section for focal orbit 



In Fig. [52] we can appreciate a pair of unstable periodic orbits { 1 }. Also included 
are a stable -j^} trajectory, a neutrally stable periodic orbit { 3 Q 6 }, as well as three 
quasi-periodic orbits. 

For the second example we start with a pair of three-bounce periodic unstable 
orbits. But now we have also included a Reuleaux type triangle (see Fig.l23|). which 
will serve as a common forbidden inner region. 

For the three neutrally stable periodic orbits in Fig. [24] as well as for the three 
quasi-periodic orbits we observe that they are confined to the region shown on the 
right hand side of Fig. [221 

To illustrate the general situation several pairs of consecutive unstable periodic 
orbits with a single quasi-periodic orbit between them are shown in Fig. 1251 

To see how all the orbits related to the {2} orbit are organized in configuration 
space see Fig. [26] 



9. Concluding remarks 

EKELAND [9]: "there is no mistaking an intcgrablc system 

with a nonintcgrable one." 

From the Poincare surface of section one can get some pretty clear idea of the 
general behaviour a dynamical system. In it we can identify 

• periodic trajectories which appear as a finite collection of points, 

• quasiperiodic trajectories which appear as one-dimensional curves 
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Figure 20. A librational, a { g 5 }, a { x 3 8 }, and a { H } orbit. 



• and finally irregular or chaotic trajectories which appear as a scatter of 
points limited to a finite area. 

In the Poincare surface of section for a typical billiard we are able to find all of these 
different types of trajectories. A Poincare surface of section consisting exclusively of 
closed curves and periodic points suggests completely regular (integrable) behavior. 

The extreme cases of completely regular and of fully chaotic behavior are both 
a very uncommon occurrence. Most often one finds that the Poincare SOS is filled 
with a mixture of both regular and irregular orbits (see Korsch and Zimmer [21] ). 

So the question is: exactly, how exceptional is the hexagonal string billiard? 

The following are the most relevant properties of this billiard: 

(1) The boundary of our billiard is twice continuously differentiable and has 
strictly positive curvature. 

(2) Our billiard is not a simple deformation or perturbation of the circular 
billiard. In addition to the planar region bounded by six elliptical arcs we 
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FIGURE 21. A self- retracing librational, a { g 8 } , a { ^ } , and a 
{ 3 9 6 } orbit. 



also have a set of very special points: the foci, whose influence cannot be 
ignored. 

(3) Our billiard exhibits three types of trajectories: focal, outer, and inner (see 
Theorem 14. 1|) . 

(4) It has a special caustic and we can associate a forbidden region with any 
trajectory. 

(5) The Poincare surface of section is filled solely with periodic points and 
invariant curves, while chaotic area-filling orbits are totally absent. 

(6) In the Poincare surface of section there are no island chains. 

Finally let us briefly mention a long-standing open problem, very often attributed 
to G. D. Birkhoff, although originally stated by Poritsky [30](see also Gutkin |15|). 
the so-called Birkhoff-Poritsky conjecture. It basically states that among all bil- 
liards inside smooth closed convex curves, only billiards in ellipses are integrable. 



NUMERICAL EXPLORATION OF A HEXAGONAL STRING BILLIARD 



2:i 




FIGURE 22. Trajectories related to the four-bounce periodic orbit 




Figure 23. Unstable three bounce orbits and Reuleaux triangle 



This means that aside from the energy, there must exist a second conserved quan- 
tity. Finding such a constant of motion is usually a very difficult task. However in 
the case of the elliptic billiard it is well-known: it is simply the product of the two 
angular momenta about the two foci (see Zhang et al. [38], Berry [4], Korsch and 
Zimmer |21j). 

So, we want to conclude with a question: 
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FIGURE 24. Trajectories related to the three-bounce periodic orbit 

from all the evidence gathered for the hexagonal string billiard, can we still be 
certain that the elliptic billiard is the only integrable convex billiard? 
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Appendix A. 

Proof of Theorem \2.1\ Consider now the ellipses £\ with foci Fi, F3 and £ n with 
F n , F2 and where the sum of the distances from the two foci is given by I — 2(n — 2) = 
2 + 2d. The length of the portion of the string that is in contact with the n-gon is 
2(n-2). 

The ellipses £\ and £ n are of course very smooth closed curves and we need 
merely verify that the equation for £ 1 has the same function value and first and 
second derivatives as the equation for £ n at x — thus ensuring C 2 continuity (see 
Fig.©. 

Letting P(x,y) be an arbitrary point, then for the ellipse £\\ 

\F X P\ + \F 3 P\ = 2 - — 
cos a 
or 



J[x + l) 2 + y 2 + y/(x- l + 2cosa) 2 + (y + 2sina) 2 = 2 

cos a 

Similarly, for the ellipse £ n : 

\F n P\ + \F 2 P\ = 2 - — 



cos a 
or 



a/0 + 1 - 2 cos a) 2 + (y + 2 sin a) 2 + \/(x- l) 2 + y 2 = 2 



cos a 

By bringing all the terms from the right-hand side to the left-hand side, we get 
the equations 

ei(x, y) = and e n (x,y) = 
Notice that since 

e n (x,y) = e x {-x,y) 
the values of the functions and of the second derivatives for the equations of both 
ellipses have to coincide when x = 0. 

A simple, but tedious process of straightforward differentiation and substitution 
shows that at the common intersection point Gi = (0, — ™12l) the first derivative is 
and the second derivative is (cos a— 1) cos a sin a/{— 1+2 cos a) for both equations. 
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□ 

Appendix B. 

Proof of Theorem \4-l\ 

i) There are two extreme supporting lines at each vertex and any ray that lies 
in the angle formed by the them is also a supporting line of K. 

So to start, without loss of generality, let us consider a billiard segment 
that lies in the angle formed by the two extreme supporting lines at vertex 
F 2 . It will have necessarily one extreme point on arc 24 and the other on 

arc 62. If the extreme point on arc 24 is chosen as the initial point Pq then 
we get a ray passing through the vertex F 2 which will then intersect the 

arc 62 at some point Pi. 

Because of the reflective property of ellipses a ray leaving one focus, in 

this case F 2 , will reflect off the ellipse corresponding to the arc 62 and pass 
through the second focus Fq. 

But then we get again a ray (through Pi and Fq) that lies in the angle 
formed by the two extreme supporting lines for the vertex Fq so it has to 
be on a supporting line itself. 

That just means that the billiard trajectory remains on a supporting 
line, 
ii) Let us now consider the angle formed by the two supporting lines of K 

through a point Pq on arc 24. One of them goes through F 2 and the other 
through P4. Because of this they are also supporting lines of the closed 
segment P2P4, which is the segment between the two foci of the ellipse 

corresponding to arc 24. Now we can consider two different types of rays 
issuing from P$: those strictly inside and those strictly outside the convex 
angle P2P0P4. The first ones will intersect the interiors of both K and of 
the segment F 2 F4(see^9\) , whereas the second ones do not intersect either 
K or the closed segment P2P4. 

In the first case we get a segment P0P1 where Pi is some point belonging 
to the arc j,j + 2. Let us consider the reverse trajectory, the one starting 
at point Pi and ending at point Pq. It, of course, intersects the interior 
of K. But then it has to intersect the interior of the segment FjFj +2 . So 
then the original trajectory intersects this segment and it has to do so again 
after reflection with the boundary. This implies that the reflected segment 
of the trajectory then also intersects the interior of K . 
iii) We leave this case to the reader. 

□ 

Appendix C. 
Proof of Theorem \4-£\ Because of the symmetry present in our billiard, we can 



restrict attention to a typical arc, say 24. 

Let P(x,y) be a point on this arc (See Fig. [3J. Then 



-12 

cos 7 = 



lab 
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where 



so 



= ^(.T-l) 2 + (y + V3) 2 , b = ^{x - l) 2 + (y- Vzy 
-12 + 2(x - l) 2 + (y + V3) 2 + (y- %/3) 2 



cos 7 



2y/(x - l) 2 + (y + VZ) 2 ^{x - l) 2 + (y - v^) 2 
Substituting the equation for the ellipse 



0r-l) 2 = 6(l-|-) 



we get 



-4 a 2 



cos 7 



+ (y + V3) 2 + ( 2 /-V3) 2 



V 6 - ¥ + (?y + ^V 6 - ¥ + (y - ^ 



or after simplifying 

Since 
it is easy to see that 

so that 



cos 7 



Z/ 2 + 9 

27 -y 2 

\/3<y<\/3 



1 1 

- < cos 7 < - 
3 - ' ~ 2 



60° < 7 < 70.53° 



Appendix D. 



□ 



Proof of Theorem \4-S\ For the notation see Fig. [HI Consider the sequences 

{si}, {<fii}, {a;}. It is clear that they are bounded: 

• 2 < Si < 4 

• tt/6 < (fi < 7r/2 

• tt/3 < a; < arccos (|) (see ThmfO]). 

We will show that {(fi},{Si} are both nondecreasing sequences. 
Subtracting 

<Pi+i + <Pi+i = 2tt/3 
from 

an + ifii + (p' i+1 = it 
we get 

< cti - 7r/3 = tpi+i - ifii 

and therefore ^ < (fi+i. 
From 

s 2 = 12 + tf - 2 • 2 V% cos <# 
and 

Si + U = 6 
we get 

4-\/3 — 6 cos v?i 



V3 



COS </?i 
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Note that the function f (w) = -^= — ^^ is nondccreasing. 

J XT ' V 3 — cos <p 

So if ^ < ip i+ i then Sj = f(ifi) < /(v»+i) = Sj+i- 

Now {(/?i} and {s^} both converge because they are monotone bounded sequences. 

But then as i — > oo 

Q!j - 7T/3 = (Pi+l - ifi -)• 



and therefore 

From 
and 
we obtain 



tt/3 



12 = s 2 + t 2 — 2siii cos«i 



•*i = 6 

12 + s 2 - 6sj 



COS OLi = — - 



Si(-Q + Si) 
So if i — > oo then s; — > s* and 



1 

2 
which yields s* = 2,4. 
Using the law of sines we get 


12 + s* 2 -6s* 

s*(-6 + s*) 




sin cti sin ^ 



12 



Vo/2 sin (/?* 



if i — > oo then 

that is 
and so finally 
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